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We investigate the effects of quenched bond randomness on the critical properties of the two- 
dimensional ferromagnetic Ising model embedded in a triangular lattice. The system is studied in 
both the pure and disordered versions by the same efficient two-stage Wang-Landau method. In 
the first part of our study we present the finite-size scaling behavior of the pure model, for which 
we calculate the critical amplitude of the specific heat's logarithmic expansion. For the disordered 
system, the numerical data and the relevant detailed finite-size scaling analysis along the lines of 
the two well-known scenarios - logarithmic corrections versus weak universality - strongly support 
the field-theoretically predicted scenario of logarithmic corrections. A particular interest is paid to 
the sample-to-sample fluctuations of the random model and their scaling behavior that are used as 
a successful alternative approach to criticality. 
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I. INTRODUCTION 

Understanding the role played by impurities on the 
nature of phase transitions is of great importance, both 
from experimental and theoretical perspectives. First- 
order phase transitions are known to be dramatically 
softened under the presence of quenched randomness [l|- 
[To| . while continuous transitions may have their expo- 
nents altered under random fields or random bonds [111 . 
[T^ |. There are some very useful phenomenological ar- 
guments and some, perturbative in nature, theoreti- 
cal results, pertaining to the occurrence and nature of 
phase transitions under the presence of quenched ran- 
domness 0, H, [H, [ll] . Historically, the most celebrated 
criterion is that suggested by Harris [111 ]. This criterion 
relates directly the persistence, under random bonds, of 
the non random behavior to the specific heat exponent 
a p of the pure system. According to this criterion, if 
a p > 0, then disorder will be relevant, i.e., under the ef- 
fect of the disorder, the system will reach a new critical 
behavior. Otherwise, if a p < 0, disorder is irrelevant and 
the critical behavior will not change. 

Pure systems with a zero specific heat exponent (a p = 
0) are marginal cases of the Harris criterion and their 
study, upon the introduction of disorder, has been of 
particular interest [l5|. The paradigmatic model of the 
marginal case is, of course, the general random 2d Ising 
model (random-site, random-bond, and bond-diluted) 
and this model has been extensively investigated and de- 
bated [see Ref. [16| and references therein]. Several re- 
cent studies, both analytical (renormalization group and 
conformal field theories) and numerical (mainly Monte 
Carlo (MC) simulations) devoted to this model, have 
provided very strong evidence in favor of the so-called 
logarithmic corrections's scenario [l7l - [2l| . According to 
this, the effect of infinitesimal disorder gives rise to a 
marginal irrelevance of randomness and besides logarith- 
mic corrections, the critical exponents maintain their 2d 



Ising values. In particular, the specific heat is expected 
to slowly diverge with a double-logarithmic dependence 
of the form C oc ln(lnt), where t = \T - T c \/T c is the 
reduced critical temperature [l7l - [2l| . Here, we should 
mention that there is not full agreement in the literature 
and a different scenario, the so-called weak universality 
scenario [22T4251 ] . predicts that critical quantities, such 
as the magnetization, zero-field susceptibility, and cor- 
relation length display power-law singularities, with the 
corresponding exponents /?, 7, and v changing continu- 
ously with the disorder strength; however this variation 
is such that the ratios (3/v and j/v remain constant at 
the pure system's value. The specific heat of the disor- 
dered system is, in this case, expected to saturate, with a 
corresponding correlation length's exponent v > 2/d [T3 ] . 

In general, a unitary and rigorous physical description 
of critical phenomena in disordered systems still lacks 
and certainly, lacking such a description, the study of 
further models for which there is a general agreement 
in the behavior of the corresponding pure cases is very 
important. In this sense, the triangular Ising ferromag- 
net is a further suitable candidate for testing the above 
predictions that has not been previously investigated in 
the literature. Thus, our investigation will be related to 
the extensive relevant literature concerning the critical 
properties of the disordered 2d square Ising model [lB - 
123. 12q - t53 ] . In particular, our discussion will focus on the 
main point of dispute over the last two decades, concern- 
ing the two well-known conflicting scenarios mentioned 
above and we will provide additional new evidence in fa- 
vor of the well-established scenario of strong universality. 
We should note here that, the theoretically predicted sce- 
nario of strong universality has been confirmed by several 
MC studies on the square lattice starting from the early 
90 's to nowadays [see Ref. ffB] for a detailed historical 
review] . 

As mentioned above, it is always important to consider 
further models for which the critical properties of the 
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corresponding pure versions are exactly known. This is 
also the case for the present model under consideration, 
namely the triangular Ising model, called hereafter as the 
TrIM, defined as usual by the Hamiltonian 

H = -jJ2^j, (1) 

<ij> 

where the spin variables s, take on the values — 1,+1, 
< ij > indicates summation over all nearest-neighbor 
pairs of sites, and J > is the ferromagnetic exchange 
interaction. The TrIM belongs to the same universality 
class with the corresponding square Ising model (SqlM) , 
sharing the same values of critical exponents and a log- 
arithmic behavior of the specific heat jfi(| [6l| . Addition- 
ally, the critical temperature of the model and also the 
critical amplitude A of Ferdinand and Fisher's [IH spe- 
cific heat's logarithmic expansion [see also the discussion 
in Sec. Mil and Eq. (01 are exactly known from the early 
work of Houtappel ]63] to be T c = 4/ln3 = 3.6409- ■ • 
and A = 0.499069 • • • , respectively. Nevertheless, it ap- 
pears that for the TrIM a verification of the finite-size 
scaling (FSS) properties of the model using high qual- 
ity data from MC simulation is still lacking. Thus, the 
first part of this work is devoted to the investigation of 
the FSS behavior of the model, especially the estimation 
of the amplitudes and other relevant coefficients in the 
specific heat's logarithmic expansion and also to the es- 
timation of the critical exponents. In this sense, the aim 
of this first part is twofold: First, to provide the first 
detailed FSS analysis of the pure model and, second, to 
present a concrete reliability test of the proposed numer- 
ical scheme. 

Our main focus, on the other hand, is the case with 
bond disorder given by the bimodal distribution 

P(Jn) = \ Wa - Ji) + s ( J v - ; (2) 

= I ; Ji > J 2 > ; r = — , 

A J\ 

so that r reflects the strength of the bond randomness 
and we fix 2ks /(J1+J2) = 1 to set the temperature scale. 
The value of the disorder strength considered throughout 
this paper is r = 1/3. The resulting quenched disor- 
dered (random-bond) version of the Hamiltonian defined 
in Eq. flTJ reads now as 

H =~Y1 J ij s *Sj ( 3 ) 

<ij> 

and will be referred in the sequel as the random-bond 
triangular Ising model (RBTrIM). The corresponding 
random-bond SqlM will be denoted hereafter respectively 
as RBSqIM. The model on the square lattice has the ad- 
vantage that the critical temperature is exactly known 
as a function of the disorder strength r by duality re- 
lations (64[. For the RBTrIM there exist only several 
approximations for the critical frontier of the site- and 



bond-diluted cases, obtained via renormalization-group 
techniques 65] and, to our knowledge, a study of the 
critical behavior of the model is lacking. 

The rest of the paper is laid out as follows: In Sec. |TT] 
we present the necessary simulation details of our numer- 
ical scheme. In Sec. IIIII we discuss the FSS behavior of 
the pure model, testing with our high accuracy numeri- 
cal data the exact expansion of the critical specific heat. 
Then, in Sec. [IV] we present a detailed FSS analysis for 
the random version of the model, including - apart from 
the classical FSS techniques - concepts from the scaling 
theory of disordered systems. Our results and the rele- 
vant discussion clearly favors the scenario of strong uni- 
versality in marginal disordered systems. Finally, Sec. IVl 
summarizes our conclusions. 



II. SIMULATION DETAILS 

Resorting to large scale MC simulations is often nec- 
essary j66j, especially for the study of the critical be- 
havior of disordered systems. It is also well known (67j 
that for such complex systems traditional methods be- 
come very inefficient and that in the last few years sev- 
eral sophisticated algorithms, some of them are based 
on entropic iterative schemes, have been proven to be 
very effective. The present numerical study of the RB- 
TrIM will be carried out by applying our recent and ef- 
ficient entropic scheme 59, 68, 69]. In this approach we 
follow a two-stage strategy of a restricted entropic sam- 
pling, which is described in our study of random-bond 
Ising models (RBIM) in 2d [59| and is very similar to 
the one applied also in our numerical approach to the 
3d random-field Ising model (RFIM) [69|. In these pa- 
pers, we have presented in detail the various sophisti- 
cated routes used for the restriction of the energy sub- 
space and the implementation of the Wang-Landau (WL) 
algorithm [70j. Further details and an up to date imple- 
mentation of this approach, especially for the study of 
disordered systems, is provided in our recent paper on 
the universality aspects of the pure and random-bond 2d 
Blume-Capel model [71] . 

We do not wish to reproduce here the details of our 
two-stage implementation and the practice followed in 
our scheme for improving accuracy by repeating the sim- 
ulations. However, we should like to include a brief dis- 
cussion on the approximate nature of the WL method. 
The usual WL recursion proceeds by modifying the den- 
sity of states G(E) according to the rule G(E) —> fG(E) 
and initially one chooses G(E) = 1 and / = /o = e. Once 
the accumulative energy histogram is sufficiently flat, the 
modification factor / is redefined as: fj + i = y/f], with 
j = 0,1,... and the energy histogram reset to zero until / 
is very close to unity (i.e. / = e 10 ~ 8 « 1.000 000 01). As 
has been reported by many authors in the study of sev- 
eral models, once / is close enough to unity, systematic 
deviations become negligible. However, the WL recur- 
sion violates the detailed balance from the early stages of 
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FIG. 1: Disorder distribution of the susceptibility maxima of 
two lattices with linear sizes L — 120 and L — 200 of the 
RBTrlM. The running averages over the samples is shown by 
the thick solid lines. 



the process and care is necessary in setting up a proper 
protocol of the recursion. In spite of the fact that the 
WL method has produced very accurate results in sev- 
eral models, it is fair to say that there is not a safe way to 
access possible systematic deviations in the general case. 
This has been pointed out and critiqued in a recent re- 
view by Janke |72j |. For the 2d Ising model (where exact 
results are available for judgement) the WL method has 
been shown to converge very rapidly. Furthermore, from 
our experience and especially from our recent study of 
the 2d RBSqIM [5{|, for which the exact phase diagram 
is known by duality relations, our high-level WL imple- 
mentation has produced excellent results, enabling us to 
discriminate between competing theoretical predictions 
on that model. Since the RBTrlM is expected to have a 
similar "entropy structure" to the corresponding square 
model, we anticipate and it will be verified in the se- 
quel, that, our WL scheme produces sufficiently accurate 
estimates enabling us, also in this case, to distinguish 
between competing theoretical predictions, as we have 
already done in the corresponding model on the square 
lattice. 

Using this scheme we performed extensive simulations 
for several lattice sizes in the range L = 20 — 200, over 
large ensembles {1, • • • , q, ■ ■ ■ , Q} of random realizations 
(Q = 500). Let us note here that for the pure model 
we simulated for each lattice size, 200 independent runs 
(WL random walks). It is well known that, extensive 
disorder averaging is necessary when studying random 
systems, where usually broad distributions are expected 
leading to a strong violation of self-averaging 
measure from the scaling theory of disordered systems, 
whose limiting behavior is directly related to the issue 
of self-averaging [73|, [ll] may be defined with the help of 
the relative variance of the sample-to-sample fluctuations 
of any relevant singular extensive thermodynamic prop- 
erty Z as follows: R z = {[Z 2 ] av - [Z} 2 av )/[Z} 2 av . Figured 



presents evidence that the above number of random re- 
alizations is sufficient in order to obtain the true average 
behavior and not a typical one. In particular, we plot 
in this figure (for lattice sizes L — 120 and L = 200) 
the disorder distribution of the susceptibility maxima x* 
and the corresponding running average, i.e. a series of 
averages of different subsets of the full data set - each 
of which is the average of the corresponding subset of a 
larger set of data points, over the samples for the simu- 
lated ensemble of Q = 500 disorder realizations. A first 
striking observation from this figure is the existence of 
very large variance of the values of x* q > indicating the ex- 
pected violation of self-averaging for this quantity. This 
figure illustrates that the simulated number of random 
realizations is sufficient in order to probe correctly the 
average behavior of the system, since already for Q w 300 
the average value of %* appears quite stable. 

Closely related to the above issue of self-averaging in 
disordered systems is the manner of averaging over the 
disorder. This non-trivial process may be performed in 
two distinct ways when identifying the finite-size anoma- 
lies, such as the peaks of the magnetic susceptibility. The 
first way corresponds to the average over disorder real- 
izations ([. . .]av) and then taking the maxima ([. . -]a V ), 
or taking the maxima in each individual realization first, 
and then taking the average ([...*]„„). In the present 
paper we have undertaken our FSS analysis using both 
ways of averaging and have found comparable results for 
the values of the critical exponents, as will be discussed 
in more detail below. Closing this brief outline, let us 
comment on the statistical errors of our numerical data. 
The statistical errors of our WL scheme on the observed 
average behavior, were found to be of relatively small 
magnitude (of the order of the symbol sizes) when com- 
pared to the relevant disorder-sampling errors (due to 
the finite number of simulated realizations). Thus, the 
error bars in most of our figures below concerning the 
average [. . .]* v and used also in the corresponding fitting 
attempts, reflect the disorder-sampling errors and have 
been estimated using groups of 50 realizations via the 
jackknife method [63|. On the other hand for the case 
[. . *]av the error bars shown reflect the sample-to-sample 
fluctuations. 



III. PURE MODEL 

In this Section we proceed to investigate the critical 
behavior of the pure TrIM defined in Eq. ([1]). Our aim 
is to observe the exact critical behavior of the model and 
also to estimate the whole set of critical exponents, pay- 
ing particular attention to the FSS behavior of the crit- 
ical specific heat. As mentioned above, the numerical 
data shown below in Figs. [H -[3] have been estimated as 
averages over 200 independent runs together with the 
corresponding error bars. 

Figure (Ha) gives the shift behavior of the pseudocrit- 
ical temperatures corresponding to the peaks of the fol- 



4 



4.2 
4.1 

4.0 
3.9 
3.8 
3.7 
3.6 



' I 1 1 

f t T = 3.6412(5) = 


1 1 1 

= T (exact) 


i 1 i 

= 3.6409 


n 111 v = 1 no4.(s>\ 


■ 


Z = 


C 




o 


z - 






A 


z = 


e<|M|>/aK 




V 


z = 


a In <M> / a K 




♦ 


z - 


3 In <M 2 > / a K - 






z = 


31n<M 4 >/gK • 


, . , 


1 




i . i 



40 



80 



120 



160 



200 



L 



ii 

CO 

1 

V 

d 



1000 p 



100 r 



10 



: 1 


t 1 — i — i — i — m 1 i ^ . 




■ n= 1: v = 0.998(4) • 




• n = 2: v = 0.999(5) 


:0>) 


* n = 4: v = 0.999(3) - 
■ ■ : 



10 



100 



FIG. 2: (a) Simultaneous fitting of the form © of 6 pseudo- 
critical temperatures defined in the text, (b) FSS of several 
powers of the logarithmic derivatives [Eq. |[5|] of the order 
parameter at the critical temperature in a log-log scale. The 
solid lines are linear fittings giving the value v = 1 for the 
correlation length's exponent. 



lowing six quantities: specific heat C, magnetic suscepti- 
bility Xj derivative of the absolute order par ameter with 
respect to inverse temperature K = 1/T [75j 



d(\M\) 
dK 



\M\H)-(\M\)(H), 



(4) 



and logarithmic derivatives of the first (n = 1), second 
(n = 2), and fourth (n = 4) powers of the order parame- 
ter with respect to inverse temperature 

d\n(M n ) _ {M n H) 



75] 



dK 



(M n ) 



-(H). 



(5) 



Fitting our data for the whole lattice range to the ex- 
pected power-law behavior 



Tl z = T C + bL- 1/u , 



(6) 



where Z stands for the different thermodynamic quan- 
tities mentioned above, we find the critical temperature 
to be T c — 3.6412(5) which is in excellent agreement 
with the exact value 3.6409 • ■ ■ . Additionally, our esti- 
mate of the critical exponent v of the correlation length 
is v = 1.004(5), also in excellent agreement with the value 
v = 1. Additional estimates for the critical exponent v 
may be obtained from the scaling behavior of the loga- 
rithmic derivatives ([5]), which scale as ~ L 1 ^ with the 
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FIG. 3: FSS of the specific heat data at the exact critical 
temperature. 

lattice size [75] . The FSS of these logarithmic derivatives 
at the critical temperature is shown in Fig. HJb) and in 
all three cases a value v w 1 is obtained consistent with 
the estimate from panel (a) of Fig. [5] and with the exact 
value v = 1. 

We know turn to the most interesting issue in the study 
of the pure model, which is the specific heat's logarith- 
mic expansion, as mentioned in the introduction. For 
the square lattice, it was shown in 1969 in the pioneer- 
ing work of Ferdinand and Fisher [62| that close to the 
critical point the specific heat obeys the following FSS 
expansion 



C L {T) = A \nL + B(T) 



1 



B 2 (T)- 



„ In L 
B l{ T) — 

where the value of the critical amplitude A is 
0.494358 ■ • • . As pointed out in Ref. [fj3], this is the same 
with the amplitude Aq in the temperature expansion of 
the specific heat close to the critical point and this was 
already known from the original paper of Onsager (76[. 
The first B coefficients are given in Ref. 62] and further 
details have been presented in Refs. [68|, I77H79I ]. In par- 
ticular, at the critical temperature the constant term B 
is 0.138149 • ■ ■ and as it is also well-known [H [zi| the 
coefficient B\ is zero. 

The universality of the above expansion has been al- 
ready pointed out and discussed by Fisher For the 
three most common 2d lattices, i.e. the square, plane 
triangular, and honeycomb, a unified approach has been 
presented by Wu et al. [79[ , from which one can find also 
for the plane triangular lattice the first B-coefficients and 
the amplitude Aq at the critical temperature. As shown 
by Wu et al. [79( for all the three lattices, the coeffi- 
cient B\ is zero at the critical temperature and as in the 
square lattice, also for the triangular lattice, the crit- 
ical amplitude of the expansion is Aq = 0.499069 • • • , 
equal to the amplitude of the temperature expansion 
close to the critical point obtained by Houtappel [63l j . 
From the work of Wu et al. [7^] we find, using their 
closed form expressions, for the plane triangular lat- 
tice (with periodic boundary conditions and aspect ratio 
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FIG. 4: Magnetic exponent ratios of the pure TrIM: FSS in a 
log-log scale of (a) the order parameter and (b) the magnetic 
susceptibility at the critical temperature. In both cases linear 
fittings are applied. 



R = 1) the following values B = B(T C ) = 0.14185- •• 
and B 2 = B 2 (T C ) = -0.15003 

It is of interest at this point to examine the compat- 
ibility of our numerical data with the above expansion 
for the case of the critical specific heat data, i.e. the 
data of the specific heat at the exact critical tempera- 
ture of the TrIM. In Fig. [3] we consider only the first two 
terms in the above expansion (i.e. we set also B 2 = 0) 
and pay attention in estimating the critical amplitude 
Aq and the constant ^-contribution. From the results 
of the fitting in the total lattice range L = 20 — 200, as 
also shown in the figure, we observe that the estimated 
value for the critical amplitude Ao = 0.4990(9) is very 
close to the exact value Aq = 0.499069 • • • [63|]. For the 
first ^-coefficient we find the value B = 0.138(5) which 
is in good agreement with the exact result 0.14185 
Let us note here that if we try to fit the data of Fig. 02 
including also the coefficient B 2 of the expansion (7]| we 
get the estimates 0.4952(59), 0.153(16), and -0.20(12) 
for the critical amplitude Aq and the coefficients B and 
B 2 , respectively. However, if we fix the values of Ao and 
B to their exact ones, the estimate for B 2 we get from 
the fitting is —0.158(11), in excellent agreement with the 
exact value, whereas if we fix the value of B and B 2 to 
their exact ones we get the estimate 0.4990(3) for the 



critical amplitude Aq. 

Finally, Figs. |4}a) and (b) present our estimations for 
the magnetic exponent ratios j3/v and j/u. For the esti- 
mation of (5/v we use the values of the order parameter 
at the exact critical temperature. As shown in panel 
(a), in a log- log scale, the linear fitting provides the es- 
timate P/v = 0.1252(4). In panel (b) we show the FSS 
of the critical susceptibility, also in a log-log scale. The 
straight line is a linear fitting for L > 20 giving the esti- 
mate 7/1/ = 1.754(5). Thus, our results presented in this 
Section for the pure 2d TrIM model are in excellent agree- 
ment with the exact results and also with the expected 
logarithmic expansion of the specific heat [HI, [79[ . This 
consists a very strong accuracy test of the proposed two- 
stage WL entropic sampling in restricted energy spaces. 



IV. RANDOM MODEL 

We now present our numerical results for the random- 
bond version of the triangular Ising model for disorder 
strength r — 1/3. From simple universality- type theoret- 
ical arguments, this system is also expected to undergo 
a second-order transition between the ferromagnetic and 
paramagnetic phases and in particular it should be also 
expected that this transition will be in the same univer- 
sality class as the RBSqIM. 

We start our analysis of the RBTrIM by presenting 
the general shift behavior of various pseudocritical tem- 
peratures of the model. Figure 03a) illustrates the shift 
behavior of 7 pseudocritical temperatures defined as the 
temperature where the corresponding average thermo- 
dynamic property attains its maximum. The first 6 
are as those defined in Sec. [TT] for the corresponding 
pure model. The last pseudocritical temperature is a 
newly introduced temperature, defined as the tempera- 
ture where the ratio R\ x »] av defined in Sec. lU becomes 
maximum. The solid lines show an excellent simultane- 
ous power-law fitting attempt of the form ([6]) giving the 
value T c = 3.4642(52) for the critical temperature of the 
random model and a value v = 0.997(6) for the critical 
exponent v of the correlation length. The fitting shown 
in Fig. 03a) has been performed for all lattice sizes and 
it was very stable when shifting the L-range to larger 
values. Note also that a simple fitting attempt using 
only the newly defined pseudocritical temperature de- 
fined with the help of the sample-to-sample fluctuations 
gives a value T c = 3.4677(47) for the critical tempera- 
ture and a value v = 0.994(8) for the correlation length's 
exponent. These overall estimates for the exponent v 
consist a strong indication that the RBTrIM shares the 
same value of v as the pure version, thus reinforcing the 
scenario of logarithmic corrections (strong universality). 

Figure OJb) illustrates again the shift behavior of the 6 
pseudocritical temperatures of panel (a), estimated now 
via the second way of averaging discussed in Sec. |H i.e. 
by taking the average over the individual pseudocriti- 
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FIG. 5: (a) - (b) Shift behavior of several pseudocritical 
temperatures defined in the text. The error bars in panel 
(b) reflect the sample-to-sample fluctuations, (c) FSS of the 
sample-to-sample fluctuations of the pseudocritical tempera- 
ture of the magnetic susceptibility shown in panel (b). 



cal temperatures. The error bars shown in this panel 
reflect the sample-to-sample fluctuations of the pseudo- 
critical temperatures. Again, the results obtained from 
the simultaneous fitting attempt T c = 3.4682(59) and 
v = 0.996(8), as also shown in the figure, agree excel- 
lently with the estimates of panel (a), providing further 
evidence in favor of the accuracy of our numerical scheme 
and the strong universality hypothesis. Noteworthy that, 
if we fix the exponent v to the exact value v = 1 we get 
the most accurate estimates for the critical temperature 
to be 3.4663(16) and 3.4669(19) from the corresponding 
fittings of panels (a) and (b) of Fig. [5j 

Using now the above sample-to-sample fluctuations of 



^ 1000 



A 
V 



100 



10 



■ n = 1: v = 


1.005(7) 


(«■): 


• n = 2:v = 


1.002(5) 




r A n = 4: v = 
■ 


1.001(6) 








L =60- 


■ ■ ■ 





min 



10 



100 



L 



b 1.00 




0.96 



0.92 



0.06 



FIG. 6: (a) FSS of the logarithmic derivatives of the order 
parameter in a log-log scale. The solid lines are simple linear 
fittings for the larger lattice sizes (L m in = 60). (b) Values 
of effective exponents v e ff obtained from the data of panel 
(a) from several fitting attempts in the range (L m in — L max )- 
The solid line marks the proposed estimate v = 1. 



the pseudocritical temperatures and the theory of FSS in 
disordered systems as introduced by Aharony and Har- 
ris p73| and Wiseman and Domany 74], one may fur- 
ther examine the nature of the fixed point that controls 
the critical behavior of the disordered system. According 
to the theoretical predictions [73|, [74| , the pseudocritical 
temperatures T% of the disordered system are distributed 
with a width <5[T^] at ,, that scales with the system size as 



6([T* z ] av ) ~ L- 



(8) 



where n = d/2 or n = l/v r , depending on whether the 
disordered system is controlled by the pure or the ran- 
dom fixed point, respectively. In panel (c) of Fig. [S] we 
plot these sample-to-sample fluctuations of the pseudo- 
critical temperature of the magnetic susceptibility. The 
solid line shows a very good power-law fitting giving the 
value 1.001(3) for the exponent n of the theory, which 
is in agreement with the case n = d/2, indicating that 
the random model is controlled by the pure fixed point. 
We should note here that analogous results to those dis- 
cussed here in panel (c) for the case of the site-diluted 
Ising model on the square lattice have been presented by 
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FIG. 7: FSS of the specific heat maxima obtained via the 
two distinct ways of disorder averaging. The solid and dotted 
lines are double-logarithmic fittings of the form © for lattice 
sizes in the range L — 60 — 200. The inset shows the data 
for the case [C*]ai> as a function of the double logarithm of 
L. The solid line is an excellent linear fitting. 



Tomita and Okabc, using the probability-changing clus- 
ter algorithm [48| . 

As in the pure case, the second alternative estimation 
of v is carried out by analyzing the divergency of the log- 
arithmic derivatives of the order parameter. In Fig. [SJa) 
we illustrate in a double-logarithmic scale the size depen- 
dence of the first- (filled squares), second- (filled circles), 
and fourth-order (filled triangles) logarithmic derivatives 
(averaged over the individual maxima). The solid lines 
show linear fittings for the sizes L > 60. In all cases 
a value v = 1 is obtained for the critical exponent u, 
providing further evidence to the strong universality sce- 
nario emerged from Fig. [3] Figure [B^b) illustrates our 
method to evaluate and discuss the stability of the esti- 
mation for the exponent v from the scaling behavior of 
the logarithmic derivatives of panel (a). It shows val- 
ues of effective exponents (v e ff) determined by impos- 
ing a lower cutoff (L m i n ) and applying simultaneous fit- 
tings in windows (Lmin — L max ), where as for the pure 
case, L m ax = 200 and L mm = 20,40,60,80, and 100 
as a function of 1/L m i n . The effective estimates show 
a finite-size effect for small values of the lower cutoff, 
whereas and for L > 60 a clear trend towards the value 
v = 1 of the Ising universality class is obtained. Let us 
note here that the same picture emerged from the FSS 
of the disorder-averaged logarithmic derivatives of the 
form [d In (M n ) / dK]* av that corresponds to the first way 
of averaging, but is omitted here for brevity. We should 
note here that a similar cross-over behavior in the esti- 
mates of the critical exponent v has been observed in the 
case of the 2d site-diluted SqlM by Ballesteros et al. [39[ 
and has been explained as logarithmic corrections. Thus, 
summarizing our estimates for the critical exponent v 1 we 
feel that it is clear that it maintains the value v = 1 of 
the pure case, indicating again the validity of the strong 
universality scenario. 



We continue the presentation of our results by showing 
in Fig. [7] the FSS of the specific heat maxima averaged 
over disorder: [C]* v (up filled triangles) and [C*] u (down 
open triangles). Using these data for the larger sizes 
L > 60, we tried to observe the quality of the fittings, 
assuming a double-logarithmic divergence of the form 



[C]:„; [C* 



Ci+C 2 ln(lnL), 



or a simple power law 



[C]* av - 1 [C*] av ^C 00 +C 3 L a ^. 



(9) 



(10) 



Although it is rather difficult to numerically distinguish 
between the above scenarios, our detailed fitting at- 
tempts indicated that the double-logarithmic scenario 
[Eq. applies better to the numerical data and this 
is generally true for both [C]* av and [C*] al , data. 

In fact, the double-logarithmic fitting is shown in 
the main panel, whereas in the corresponding inset of 
Fig. [7] the data of [C*] at) are plotted as a function of 
ln(lnL). The solid line shown is an excellent linear fit 
for L > 60. Let us now give some details on the qual- 
ity of the applied fittings. We used the following sets 
of data points (L min - L max ), with L max = 200 and 
L mm = 20, 40, 60, 80, and 100. The quality of the fittings 
indicated a very good trend for the values of x 2 /DoF for 
the double logarithmic fittings ^ in the range: 0.2 — 0.7 
and for both sets of data points. However, a strong relia- 
bility test in favor of the logarithmic corrections scenario 
is provided by the stability of the coefficient C2, for both 
[C]* av (C2 « 1.43(5)) and [C*] av {C 2 « 1.48(4)) data. On 
the other hand, the estimated values of the exponent a/v 
of the power law (JTU]), for both [C}* av and [C*] at ,, fluctu- 
ate in the range [—0.12(9), —0.05(6)] (as L min increases) 
with the fitting procedure becoming rather unstable as 
we move to larger values of L m i n . The conclusion is that 
our numerical data are more properly described by the 
double logarithmic form ([9]), in agreement with the MC 
findings of Selke et al. (43] and Ballesteros et al. 
for the site-diluted SqlM and also with those of Wang 
et al. [28j for the strong disorder regime (r = 1/4 and 
r = 1/10) of the RBSqIM. 

In Fig. [8] we provide estimates for the magnetic expo- 
nent ratios j3/v and 7/1/ of the RBTrlM. In panel (a) we 
plot the average magnetization at the estimated critical 
temperature, as a function of the lattice size L in a log- 
log scale. The solid line is a linear fitting for L > 20 
giving within error bars the value of the pure model, i.e. 
fi/v = 0.1253(5) w 0.125. Additional estimate for the 
ratio (3/v can be obtained from the FSS of the derivative 
of the absolute order parameter with respect to inverse 
temperature defined in Eq. which is expected to scale 
as L^ 1 "' 3 )/ 1 ' with the system size [7a]. Thus, in panel (b) 
of Fig. [5] we plot the data for d{\M\) jdK averaged over 
disorder as a function of L, also in a double-logarithmic 
scale. The solid line is a linear fitting for the larger lat- 
tice sizes L > 60, which combined with the value v = 1, 
gives an estimate of 0.1247(4) for the ratio /3 jv. Finally, 
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in panel (c) we present the FSS of the maxima of the 
average magnetic susceptibility [x]a V (filled squares) and 
also the average of the individual maxima [x*]av (open 
triangles). The solid and dotted lines present linear fit- 
tings using the total lattice range spectrum, giving the 
estimates 1.751(5) and 1.756(9) for the ratio j/v in very 
good agreement with the expected value 1.75 of the pure 
system. For the average [x]av ^ ne error bars indicate the 
statistical errors due to the finite number of the realiza- 
tions, as discussed in Sec. Q] For the average [x*]av the 
errors bars shown reflect now the relatively large samplc- 
to-samplc fluctuations. 

Finally, using the latter sample-to-sample fluctuations, 
we construct the ratio R[ x *] av an d plot it as a function of 
the inverse linear size, as shown in the inset of Fig. [Etc). 
Clearly, for the present model the limiting value of R[x*] av 
is non-zero, indicating, as expected also for marginal dis- 
ordered systems [73], a strong violation of self-averaging. 



V. CONCLUSIONS 

The effects induced by the presence of quenched bond 
randomness on the critical behavior of the 2d Ising spin 
model embedded in the triangular lattice have been in- 
vestigated by an efficient entropic scheme based on the 
Wang-Landau algorithm. In the first part of our study 
we presented the finite-size scaling behavior of the pure 
model, for which we calculated with high accuracy the 
critical exponents and the coefficients of the specific 
heat's logarithmic expansion at the critical point. Our 
results are in full agreement with the exact expansion 



presented by Wu et al. [79(. 

In the second part of our study we investigated the 
critical properties of the disordered system. The pre- 
sented detailed finite-size scaling analysis along the lines 
of the two existing scenarios - strong versus weak univer- 
sality - strongly supports the scenario of strong universal- 
ity. Thus, our results are in agreement with the behavior 
predicted originally on theoretical basis many years ago 
by Dotsenko and Dotsenko 11711 ■ Jug 181, Shalaev [19j . 
Shankar [20|, and Ludwig |2l| and verified by simula- 
tions in recent years for the square random Ising model 
by several authors [2§El El 13 M, HI . Particular 
interest was paid to the sample-to-sample fluctuations of 
the random model and their scaling behavior that were 
used as a successful alternative approach to estimate the 
critical temperature and the correlation length's expo- 
nent. Closing, we would like to note that another inter- 
esting candidate, that has not been studied before in the 
triangular lattice, is the 3-state Potts model, which, in its 
pure version, has a positive specific-heat exponent. Dis- 
order would be relevant in this case and could provide a 
further complementary study of the present work, analo- 
gous to the early transfer-matrix calculations of Derrida 
et al. [8(| of the random-bond model on the square lat- 
tice. 
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FIG. 8: (a) FSS of the average order parameter at the crit- 
ical temperature, (b) FSS of the disorder-averaged inverse- 
temperature derivative of the absolute order parameter, (c) 
FSS of the disorder-averaged magnetic susceptibility. The in- 
set shows the limiting behavior of the ratio R[ x *] av - In au 
main panels (a) - (c) a double logarithmic scale is considered. 
Additionally, all fitting attempts are performed in the com- 
plete lattice range L = 20 — 200 since no deviation in the 
estimated values of the critical exponents was observed when 
shifting the value of the lower cutoff L min . 



